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Abstract 

The binding energies of deformed even-even nuclei have been analysed within the framework 
of a recently proposed microscopic-macroscopic model. We have used the semiclassical Wigner 

- Kirkwood h expansion up to fourth - order, instead of the usual Strutinsky averaging scheme, 
to compute the shells corrections in a deformed Woods - Saxon potential including the spin-orbit 
contribution. For a large set of 561 even-even nuclei with Z > 8 and N > 8, we find an rms 
deviation from the experiment of 610 keV in binding energies, comparable to the one found for the 
same set of nuclei using the FRDM of Moller and Nix (656 keV). As applications of our model, we 
explore its predictive power near the proton and neutron drip lines as well as in the superheavy 
mass region. Next, we systematically explore the fourth - order Wigner - Kirkwood corrections 
to the smooth part of the energy. It is found that the ratio of the fourth - order to the second 

- order corrections behaves in a very regular manner as a function of the asymmetry parameter 
/ = (N — Z)/A. This allows to absorb the fourth - order corrections into the second - order 
contributions to the binding energy, which enables to simplify and speed up the calculation of 
deformed nuclei. 

PACS numbers: 21.10.Dr, 21.60.-n 
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I. INTRODUCTION 



The models of nuclear masses are continuously challenged by the advances in experimental 
techniques which nowadays are extending the nuclear chart to previously unexplored regions 
of exotic isotopes and superheavy elements. The theoretical description of nuclear masses 
takes place primarily along two main approaches. On the one hand, in the microscopic 
nuclear models, the nuclear binding energy is obtained from calculations with energy density 
functionals based on effective nuclear interactions In the microscopic-macroscopic 



fmic-mac) models [2|, |4|, |5[, the nuclear binding energy is obtained as the sum of a part that 
varies smoothly with the number of nucleons plus an oscillatory correction originated by 
the quantum effects. The smooth part of the mic-mac models is obtained from a liquid- 
drop model approach, whereas the shell correction is usually evaluated by the Strutinsky 
averaging method in an external potential well. 

In our previous works , we have demonstrated that the Strutinsky average can 
be replaced by the semiclassical energy computed by means of the Wigner - Kirkwood 
(WK) h expansion of the one - body partition function js-15], in order to evaluate the shell 
corrections of a system of N neutrons and Z protons at zero temperature in an external 
potential. There are some reasons supporting this choice as we have discussed in Ref. jo]. 
On the one hand, it has been shown that the Strutinsky level density is an approximation 
to the WK level density 16|. On the other hand, the WK ^.-expansion of the density matrix 
has a variational content and it is possible to establish a variational theory based on a strict 



^-expansion |l5l. Il7j|. We shall point out that the WK expansion is also well suited to deal 
with nuclei close the drip lines. Although the WK level density exhibits a well known e^ 1 ^ 2 
divergence as e — > for a potential that vanishes at large distances, integrated moments of 
the level density, such as the energy and the accumulated level density, are well behaved in 



the e — > limit as it has been demonstrated in Ref.[15|. It has been shown that these shell 
corrections, along with a simple six parameter liquid drop formula, yield a good description 
of ground - state masses of spherical nuclei spanning the entire periodic table p]. The model 
has also been applied to calculate the binding energies of few deformed nuclei, with a good 
degree of success jj]. In the present work, we extend the work reported earlier [g] to the 
deformed nuclei and explore the predictions of the model in exotic scenarios such as drip 
line nuclei and the superheavy region. In this work, we mainly restrict our attention to the 
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even - even nuclei. 

One of the important conclusions of Ref. {f| is that in this model it is necessary to carry 
out the WK expansion up to the fourth - order in h to obtain accurate shell corrections, which 
implies that in this case one needs to work out derivatives of the single particle potentials 
(nuclear potential, Coulomb potential as well as the spin - orbit potential) up to the fourth 
- order, which is a rather cumbersome task. Therefore, this gives rise to an interesting and 
important question: can the effects of the fourth - order corrections to the binding energy 
be absorbed into the second - order ones? This question is important from theoretical as 
well as practical point of view. Theoretically, this would imply that the WK series has been 
partially re-summed, whereas from a practical point of view, it implies that it is sufficient to 
expand the one - body partition function up to second - order in h to obtain shell corrections 
with comparable accuracy. 

The absorption, if possible, would imply that there is a factor (we denote the factor by a), 
which may be a function of mass number, charge number, neutron number or combinations 
thereof, defined as 

E(h A ) 

such that 

E (h 2 ) + E (h 4 ) = aE (h 2 ) (2) 

where, E (h 2 ) and E(h A ), respectively, are second and fourth - order WK corrections to 

energy. This is an important issue discussed in the present article. 

We summarise the essential details of the semiclassical Wigner - Kirkwood expansion 

of the one - body partition function in the second section. The detailed results and their 

analysis forms the subject matter of the third section. The parameters of the macroscopic 

part of our mic-mac model, which also includes curvature correction |5l and the Wigner 
ft 

term [5|, have been obtained by minimizing the x 2 value of the energies using a selected 
set of 561 even-even deformed and spherical nuclei. The ability of this mic-mac model to 
describe nuclei in the exotic scenarios is explored in section 4. On the one hand, masses of 
very proton rich nuclei, measured recently [isj], are compared with the predictions of our 
model. On the other hand, the upper limit of the outer crust in neutron stars is studied, 
which involves nuclei near the neutron drip line. Finally, we explore the superheavy region, 



and compare the theoretical alpha decay Q values and the corresponding half lives with the 
experimental values 19] . The systematic investigation of the absorption factor a as defined 
above is contained in the fifth section. The summary and conclusions are given in the last 
section. 



II. FORMULATION 

For a system of N non interacting Fermions at zero temperature in a given external 
potential, the quantal one - body partition function is given by: 

Z(f3) = Tr(exp(-/3#)) • (3) 

The Hamiltonian of the system (H) is expressed as: 

H = ^V 2 + V(f) + V LS {7) , (4) 
2m 

with V(r) being the one-body central potential and Vls(j) the spin-orbit interaction. The 
replacement of the Hamiltonian in the above equations by the corresponding classical Hamil- 
tonian leads to the well - known Thomas - Fermi equations for particle number and total 
energy. The Wigner - Kirkwood semiclassical expansion amounts to expansion of the quan- 
tal one - body partition function in the powers of Planck's constant, h, yielding systematic 
corrections to the Thomas - Fermi energy and particle number 8H13]. 

As stated before, in this work, we carry out the WK expansion up to the fourth - order 
in h. With the spin - orbit interaction, the WK expansion of the partition function can be 
written schematically as: 

2®M = Z^(P) + 4to) ■ (5) 

where, Z^(/3) (Z^q(/3)) is the WK partition function for the central potential (spin - orbit 
part). The explicit expressions for these partition functions can be found in jf], [lo|. 

The level density gwx, the particle number N and the energy Ewk are obtained by 
appropriate Laplace inversion of the WK partition function, as follows: 

g WK (e) = C- e l Z$ K (P) , (6) 
N = (7) 



and 



WK 



XN 



I (32 



(8) 



Here, A is the chemical potential, determined to ensure the correct particle number. 

The focus of the present article being the WK energy, we present the explicit expressions 
for the WK energies alone. Following Jennings et al. 10], the energy (Eq. (JSJ)) can be 
written as: 



E WK = XN- (Eg* + Eg* + E^4 N ) - (EiP + E s n P) 



(9) 



where, E9 k N denote the contribution to the energy of the order h k arising from Laplace 
inversion C^ 1 (Z^((3)/(3 2 y On the other hand, E*jjf? are corrections to the energy of the 
order h k due to Laplace inversion C^ 1 ^Z^,(/3)//3 2 j. The explicit expression are as follows 
(see js] for further details): 
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In these expressions, V is the mean field, / is the spin - orbit form factor, k is the strength 
of spin - orbit interaction, and A is the chemical potential. 

The shell corrections, which are the difference between the quantum mechanical and 
the corresponding averaged energies, can now be obtained by subtracting Ewk from the 
quantum mechanical energy. For our calculations we choose a Woods-Saxon potential as 



mean field and a suitable Woods - Saxon form factor in the spin - orbit sector. These 
potentials are generalised for taking into account deformation effects and their corresponding 

las been obtained by folding the 
6] . In the microscopic part we have 



parameters are given in Ref. Jfj]. The Coulomb potentia 
proton density distribution with the Coulomb interaction 



also included pairing correlations using the Lipkin - Nogami scheme 2CH22I]. as described in 
details in Ref. js]. 



III. CALCULATION OF BINDING ENERGIES 

In the present work, we generalise the liquid drop formula employed in [6| by adding a 
deformation dependent curvature energy term and the Wigner term. The curvature energy 
term is found to be important in improving the agreement achieved between calculations 
and the corresponding experimental binding energies [5j. The Wigner term is expected to 
be important for light nuclei as well as to describe nuclei close to the proton drip line. 
Therefore, the modified liquid drop formula used in this work reads: 
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(15) 



where the terms respectively represent: volume energy, surface energy, curvature energy, 
Coulomb energy, correction to Coulomb energy due to surface diffuseness of charge distri- 
bution and the Wigner energy. The coefficients k v , k s , k cur} tq and C4 are free 
parameters; T z is the third component of isospin, and e is the electronic charge. 

Several parametrisations of the Wigner term are available in the literature (see, for ex- 



ample 



23)). Here, we adopt the following ansatz for the Wigner term with a cut off 



on charge and mass numbers: 

E w = wi exp 



-w 2 



N - Z 



A 



S(Z- 20) (A - 40) 



(16) 



where, W\ and u>2 are free parameters. The cut offs on charge and mass numbers have been 
introduced since it is expected that the Wigner term will make significant contributions for 
nuclei with low masses. 

The Coulomb, surface and curvature terms appearing in the liquid drop formula, as 
defined above in Eq.f lT5|) . need to be modified for the deformed shapes. In particular, the 



Coulomb term is multiplied by 



J~c — n / tz — —rdrdr 1 



32n 2 J \f — r 1 
-15 



r-f*\dS-dS' (17) 



64tt 2 

where, the symbols have their usual meanings. Notice that the integrals have been carried 
out over nuclear volume, and the lengths have been measured in units of the radius param- 
eter R Q of the nucleus with zero deformation. The transformation from six dimensional to 
four dimensional integrals has been accomplished by following the technique developed by 
Kurmanov et al. 24|. The surface term, on the other hand, is simply modified by the ratio 
of deformed to the corresponding spherical surface areas. The curvature energy term, too, 
needs to be modified to take the deformation effects into account. The modified curvature 



E cur = / (J_ + i_ ) dS (18) 



energy (E cur ) reads: 

E I 

where, E® ur is curvature energy at zero deformation; R\ and R 2 are the principal radii of 
curvature of the nuclear surface (in the units of R a ), defined by r = r s ; and dS refers to 
the area element of the nuclear surface. The surface parametrisation assumed in the present 
work is given by: 

r s = CRo(l + J2 a ^ Y ^ ■ ( 19 ) 

Here, the Y\ tfl functions are the usual spherical harmonics and the constant C is the volume 
conservation factor (the volume enclosed by the deformed surface should be equal to the 
volume enclosed by an equivalent spherical surface of radius R ): 

-1/3 

(20) 
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X,ti J 



The term Z 2 /A, which is the correction to Coulomb energy due to surface diffuseness of the 
charge distribution, does not have any explicit deformation dependence. This is because 
the distance function chosen here is such that the surface thickness is the same in all the 
directions (see discussion about this in Ref. jf]]). 

The total binding energy of a nucleus with iV neutrons, Z protons and deformation 
parameters fa, Ai and 7 is given by: 

E{N,Z,fa,fo,i) = E LDM (N,Z,(3 2 ,fa,j) + ri6E(N,Z,fo,fa n ) (21) 

8 



where, 5E represents the microscopic part of the binding energy (shell correction plus pairing 
energy). The microscopic part has been multiplied by a factor 77, which is chosen to be 0.85. 
One of the reasons for introducing such a factor is that the Coulomb potential used in the 
present work is less repulsive near r = than the corresponding value obtained by using the 
hard sphere approximation, used in the fit of proton mean field (see discussion on this point 
inRef. fl). 

The free parameters of the liquid drop formula are determined by minimising the x 2 value 



in comparison with the experimental binding energies 251 ]: 

'E(N V Z,)-E { £ pt 



x 2 

n 



1 



AE (i) 

L - 1Lj expt 



(22) 



where E(Nj,Zj) is the calculated total binding energy for the given nucleus, E^J pt is the 
corresponding experimental value 25], and AE^J pt is the uncertainty in E~J pt . In the present 



fit, for simplicity, AE^J pt is set to 1 MeV. 

To obtain these parameters we proceed as follows. We start by setting in the liquid 
drop mass formula (115j) the values obtained in our spherical calculation jf|. Explicitly, these 
values are: a v = -15.841 MeV, a s = 19.173 MeV, k v = -1.951, k s = -2.577, r = 1.187 fm 
and C4 = 1.247 MeV. Next, we choose a set of 561 even-even nuclei with Z > 8 and N > 8, 



the list of which may be found at [26j . This set comprises doubly magic, semi magic as well 
as open shell nuclei, many of which are expected to be deformed. The main task now is to 
determine the liquid drop parameters as well as the optimal deformation parameters. The 
calculation proceeds in the following steps: 

1. Assuming the previously reported |6j values of the liquid drop parameters, the binding 
energies of these nuclei are obtained by minimising on a range of 02 values (0^ is set 
to zero in this step). This gives a preliminary estimation of 02- Next, keeping this 02 
fixed, 04 is varied to obtain minimum energy. Thus, we now have preliminary values 
of both the deformation parameters. 

2. In the next step, keeping the deformation parameters fixed as obtained in the earlier 
step, the liquid drop parameters are fitted by minimising x 2 ■ 

3. With the new values of liquid drop parameters, the deformation parameters are ob- 
tained once again as described in step 1, followed by a final re-fit to the liquid drop 
parameters. 



The numerical values of the new constants of the liquid drop formula obtained through 
this minimisation procedure are: a v = -15.435 MeV, a s = 16.673 MeV, a cur = 3.161 MeV, 
k v = -1.874, k s = -2.430, k cur = (see discussion below), r = 1.219 fm, C 4 = 0.963 MeV, 
Wi = -2.762 MeV and Wi = 3.725. The values of volume, surface and Coulomb coefficients 
differ from those reported earlier {(| , primarily due to the inclusion of curvature and Wigner 
terms and the deformation effects. The curvature term, as described earlier, depends on the 
mean curvature of the nucleus, which is a function of the geometry of the nuclear surface. 
Therefore, the curvature energy, a priory, is expected to modify the surface energy term as 
well as the Z 2 /A term, which is the correction due to the surface diffuseness of the charge 
density term. The somewhat smaller value of the volume coefficient reported here, is not 
surprising. The reduction is due to the influence of the curvature term, as it has also been 
found by Pomorski and Dudek (see Table 1 of Ref. j|5]). 

It is to be noted that the coefficient of the isospin dependent term in the curvature energy 
is very difficult to determine with experimental masses. In our case the resulting statistical 
error in the corresponding parameter turns out to be more than 50% of the numerical value 
of the coefficient. Further, this term is found to weaken the strength of the isospin dependent 
term in the surface energy by a factor of 5. The isospin dependence in the curvature term, 
therefore, has been dropped from the present investigation. 

The rms deviation of the calculated binding energies with respect to the experiment 
obtained is 610 keV. The Moller - Nix calculations [27], for the same set of nuclei, yield 
a deviation of 656 keV. The explicit values of binding energies of our selected set of 561 
even-even nuclei used in the minimisation procedure can be found at |26 ] . The present 
calculation establishes that our model is indeed capable of reproducing binding energies of 
deformed nuclei as well, with excellent accuracy. The difference between the calculated and 
the corresponding evaluated 25j binding energies is presented in Fig. [TJ The corresponding 
differences obtained for the Moller - Nix calculations is presented in the same figure for 
comparison. The excellent agreement found between the calculations and experiment is 
amply clear from the figure. 

We next present and discuss the results obtained for Sr, Sn, Gd and Po isotopes as 
illustrative examples. The difference between the fitted and the corresponding experimental 
binding energies for these nuclei are plotted in Figs. [2] and [31 along with the corresponding 
differences obtained from the Moller - Nix calculations [27j. The figures reveal that the 
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FIG. 1: Difference between the calculated (fitted) and the corresponding experimental 25] binding 
energies, as a function of mass number. The dashed horizontal lines correspond to 8BE = 610 
keV. The corresponding differences obtained by using the Moller - Nix binding energies are also 
presented for comparison. 

calculated binding energies (denoted by WK) are quite close to the experimental values. The 
differences are found to vary quite smoothly as a function of mass number. Next, we present 
the two - neutron separation energies for these chains. The two - nucleon separation energies 
highlight the shell structure in an isotopic chain. Correct prediction of these separation 
energies is crucial for determination of the drip lines. The calculated and the corresponding 
experimental 25j two - neutron separation energies are plotted in Figs. H] and El The figures 
reveal that the present calculations reproduce the experimental separation energies very well 
and that the shell - gaps are also reproduced nicely. 

In addition, the systematics of deformation parameters obtained in these calculations 
turns out to be reasonable. As an illustrative example, we focus on the Sr - Zr region. It 



is well known from the systematics of experimentally measured charge radii [28[ that the 
charge radii increase dramatically by 2% for 97 Rb, 98 Sr and 100 Zr, in comparison to their 
respective lighter isotopes. This jump may be attributed to the possibility of onset of highly 
deformed shapes in the ground - state, around this neutron number (see, for example, 29|). 
Our calculations, too, reveal existence of highly deformed ground - states (with (3 2 ~ 0.3) 
around neutron number 60, in the Sr - Zr region. The values of obtained in this work for 
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FIG. 2: The difference between the calculated and the experimental 25] binding energies for Sr 
and Sn isotopes. 



2.5 
2.0 
1.5 
1.0 

H 0.5 

a 

-0.5 
-1.0 



-1.5 



Gd - Isotopes 




WK 

Q-O MN 



2.5 
2.0 
1.5 
Q 1.0 

§, 0.5 

-0.5 
-1.0 
-1.5 



138 144 150 156 162 
Mass Number 



~i — 1 — 1 — i — 1 — 1 — i — 1 — 1 — r 
Po - Isotopes 



WK 

O O MN 




_L_L 



j—L 



186 192 198 204 210 216 
Mass Number 



FIG. 3: The difference between the calculated and the experimental [25(] binding energies for Gd 
and Po isotopes. 



Kr, Sr, Zr and Mo chains are plotted in Fig. [6j The sudden change in the ground - state 
deformation around the neutron number 60 is very clear from the figure. 

Further, it is also well known that the ground - states of 72 Kr, 76 Sr and 80 Zr have very 
large (~ 0.4) deformation. This is known to be due to population in the intruder 1*79/2 
state. Thus, the ground - state of 80 Zr is a 12 particle - 12 hole state, which is manifested 
again by an extremely large stable deformation in the ground - state of 80 Zr. This has been 
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FIG. 4: The calculated and the experimental 25[ two - neutron separation energies for Sr and Sn 
isotopes. 
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FIG. 5: The calculated and the experimental [25[] two - neutron separation energies for Gd and Po 
isotopes. 



verified independently, for example, by the relativistic mean field calculation 30|], density 
dependent Hartree Fock calculation with Skyrme interaction [31], as well as by the Hartree 
Fock band mixing calculation 32|]. The deformation parameters reported in the Moller - 
Nix table 23J, too, are consistent with the discussion above. It is gratifying to note that the 
present calculations, indeed, yield (3 2 — -0.36, -0.41 and 0.44 respectively, for 72 Kr, 76 Sr and 
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FIG. 6: Deformation parameter /32 for Kr, Sr, Zr and Mo isotopes. 



3 Zr, which is in tune with the mean field as well as the mic - mac Moller - Nix calculations 



cited above. 



IV. APPLICATIONS OF THE PRESENT MODEL TO NEAR DRIP LINE NU- 
CLEI AND SUPERHEAVY NUCLEI 

We next test the ability of the present model to describe binding energies of the neutron 
rich and neutron deficient nuclei, as well as of the superheavy nuclei. To this end, we now 
present a few exploratory calculations. 



A. Proton Drip-line Nuclei in the Ge - Kr Region 

The masses of 63 Ge, 65 As, 67 Se and 71 Kr have recently been measured [18|. These nuclei 
are very proton rich, and are expected to be close to the drip - line. Notice that these 
nuclei are odd - even and even - odd. In this preliminary test of our model near proton drip 
line, we use the simple uniform filling approach for the calculation of the pairing energy. 
The calculated binding energies and one proton separation energies (S p ) for these nuclei, 
along with the corresponding experimental values 



Nix 



181 ] and those reported by Moller and 



271 ] are presented in Table 1. The binding energies as well as S p values obtained in the 



present work are found to be quite close to the experiment. This indicates that the present 
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TABLE I: The binding energies and one proton separation energies for proton rich nuclei. 'Calc.' 
(MN) represent the results obtained in the present work (by Moller and Nix [27]). The experimental 
binding energies have been obtained from mass excess values reported by Tu et al. The 
experimental S p values have been also been adopted from Ref. 18]. 





h 


Pi 


Binding Energy (MeV) 


S p (MeV) 


Calc. 


MN 


Expt. 


Calc. 


MN 


Expt. 


63 Ge 


+0.200 


-0.010 


-529.795 


-529.266 


-530.327 


2.557 


3.315 


2.210 


65 As 


+0.210 


-0.030 


-545.168 


-544.642 


-545.699 


0.633 


0.124 


-0.090 


67 Se 


+0.220 


-0.050 


-560.598 


-560.158 


-560.698 


2.379 


3.364 


1.852 


71 Kr 


-0.330 


0.010 


-592.047 


-591.219 


-591.150 


2.304 


3.093 


2.184 



model extrapolates reliably up to the proton drip lines. The nucleus 65 As is reported to be 
slightly unbound against proton emission with S p = —90 ± 85 keV Our calculation, on 
the other hand, yields a positive value of S p for 65 As, indicating a proton bound nucleus. 
However, it should be noted that the separation energies are obtained by taking differences 
of the relevant binding energies, and hence are very sensitive to the precise details of the 
same. The fact that the theoretical separation energies obtained in this work differ from the 
corresponding experimental values only by a few hundred keV's is quite remarkable. 



B. Composition of the Outer Crust of Neutron Stars 

The masses of very neutron-rich nuclei are particularly interesting for some astrophysical 
calculations. We next compute the composition of the outer crust of a neutron star as a 
further application of our present mass model. As one moves from the surface of a neutron 
star to its interior, the outer crust is the region comprising matter at densities between ~ 10 4 
g/cm 3 and ~ 10 11 g/cm 3 . Matter at those densities consists of fully-ionised, neutron-rich 
atomic nuclei that arrange themselves in the lattice sites of a Coulomb crystal embedded in 



a degenerate electron gas [33f| . The neutron excess of the nuclei in the outer crust becomes 
larger with increasing matter density until neutron drip starts taking place at a density 
about 4 x 10 11 g/cm 3 . At that point, one leaves the outer crust and enters the so-called 
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FIG. 7: Predicted composition of the outer crust of a neutron star as a function of the density. The 
upper line depicts the variation of the neutron number N, while the lower line depicts the variation 
of the proton number Z. The composition obtained by using the Moller-Nix mass formula is also 
presented for comparison. 

inner crust of the neutron star, where the atomic nuclei are immersed in an electron gas and 
a neutron gas. 



In order to compute t 



as described in Refs. 



34 



le composition of the outer crust we follow the usual formalism 



361 ] and references quoted therein. That is, we consider cold and 



electrically neutral matter which is assumed to be in thermodynamic equilibrium and in its 
absolute ground - state. We calculate the Gibbs free energy of this system by adding the 



34j-l36| and, finally, we evaluate 



contributions of the nuclear, electronic, and lattice terms 
the equilibrium composition (Z,N) at a certain pressure by minimising the obtained Gibbs 
free energy per nucleon. 

We display our predictions for the equilibrium nuclear species present in the outer crust 
in Fig. [7J We perform the calculations within the range p = 10 7 g/cm 3 to p = 3 x 10 11 
g/cm 3 . The variation of the neutron and proton numbers with increasing crustal density 
shows a structure of plateaus that are interrupted by abrupt jumps in the composition. As 



exemplified by the N = 50 plateau, the prevalence of a given nucleon number over a large 
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range of densities is related with the shell effect due to the filling of a nuclear shell. The 
N = 50 neutron plateau also is very illustrative of the fact that, with increasing density, it 
is energetically favorable for the nuclei of the crust to capture electrons from the degenerate 
electron gas. This results in increasingly neutron-rich nuclides along the neutron plateau. 
Eventually, the mismatch between the neutron and proton numbers is too large and the 
jump to the next neutron plateau takes place in an effort to reduce the penalty imposed on 
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36]. 



the system by the nuclear symmetry energy 

At low crustal densities up to about p = 7 x 10 8 g/cm 3 , our calculations sequentially 
favor 2qFb, j^Fe, and ||Ni as the equilibrium nuclides (with ||Cr occurring in a short density 
interval between |gFe and pFe). Once the jump to the N = 50 plateau ensues at a density 
p ~ 7 x 10 8 g/cm 3 , our model predicts the sequence of increasingly neutron-rich isotones 
8 gKr, f^Se, fjGe, §|jZn, and ^Ni. After the ||Ni nucleus, it is unfavorable to move further 
to 2(;Fe and at a density p ~ 1.2 x 10 11 g/cm 3 we find that the composition of the crust 
jumps to the N = 82 plateau (where our calculations predict the occurrence of the isotones 
42 4 Mo and \l 2 Zr). We display the results obtained with the Moller-Nix mass table 23[ in 
the same Fig. [7] for comparison. Though the overall pattern is quite similar to the results 
obtained with our calculated masses, the Moller-Nix mass table predicts more structure in 
the variation of the neutron and proton numbers with the crustal density, and the jump 
to the N = 50 plateau is delayed to a little higher density. This fact suggests that in the 
present mass region the shell effects due to the filling of nuclear shells and sub-shells are 
somewhat weaker in the Moller-Nix mass formula than in our model. 



C. Superheavy Nuclei 



Production and study of 
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superheavy nuclei is of current interest from both theoretical 



40] and experimental |19l . |41| aspects. With the advent of increasingly sensitive detection 
methods, it is possible to identify the superheavy elements, and measure a decay Q values 
precisely. The elements with Z = 118 have been produced so far 19) . Here, we apply our 
mic-mac model to a few recently reported superheavy nuclei jljj. In particular, we focus on 
the a decay Q values (Q a ). The binding energies of the parent as well as the daughter nuclei, 
necessary to obtain the Q a values, are obtained within our mic-mac model by minimising 
over the deformation (/?2, Pa) mesh. The binding energy of the a particle is adopted from 
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the Audi - Wapstra compilation 25] . The calculated (Calc.) as well as the experimental Q 



values ll9( are presented in Table 2. We find that the calculated Q a values are very close to 
the experiment. This is quite encouraging, since as in the case of the separation energies, 
the Q values as well are obtained by taking differences between two large quantities. 



TABLE II: The a decay Q values and half lives (7i/ 2 ) for some of the superheavy nuclei. 







Q a (MeV) 


Tl/2 


z 


A 


Calc. 


Expt. 


Calc. 


Expt. 


118 
116 
116 


294 
293 
292 


11.76 
10.59 
10.66 


11.81±0.06 
10.69±0.06 
10.80±0.07 


0.56 ms 
136 ms 
89 ms 


0.89±^ ms 
61^ ms 

18l^ 6 ms 


116 


291 


10.89 


10.89±0.07 


22 ms 


18+f ms 


115 


288 


10.49 


10.61±0.06 


129 ms 


87^° 5 ms 


115 


287 


11.38 


10.74±0.09 


0.69 ms 


32±{f ms 


114 


289 


9.91 


9.96±0.05 


2.7 s 


2 - 6 ±a7 s 


114 
114 
113 
113 
111 


288 
287 
283 
282 
280 


10.26 
10.19 
10.82 
10.99 
9.33 


10.08±0.06 
10.16±0.06 
10.26±0.09 
10.78±0.08 
9.87±0.06 


0.28 s 
0.43 s 
4.6 ms 
17 ms 
17 s 


0.80±°;^ s 
0.481™ s 
lOOt^o ms 
73t^ 4 ms 
3.61*1 s 


111 


279 


10.56 


10.52±0.16 


5.5 ms 


170±^ ms 



The a decay Q values can be related to the half lives through the Viola - Seaborg relation 



42J. In particular, following Oganessian |l9|, we adopt: 



logT 1/2 = a + cZ + d (23) 

where, Z is the charge number of the parent nucleus; Q a is the a decay Q value, and a, b, c 
and d are parameters, taken to be [xq! : a = 1.787, b = -21.40, c = -0.2549 and d = -28.42. The 
half lives obtained by using the calculated Q values are found to be in reasonable agreement 
with the experiment. At places, the calculations do deviate by an order of magnitude, but 
notice that the half lives have very large uncertainties. 
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V. SYSTEMATIC INVESTIGATION OF THE FACTOR a 



Large scale calculations using the proposed mic-mac model can be cumbersome and highly 
time consuming. Therefore, it may be very useful to look for simplifications that allow to 
speed up the calculations without loss of accuracy. To this end, we explore the possibility 
of absorbing the fourth - order correction 

E 4 = Eg N (n) + Eg?(n) + Eg N (p) + E s h °(p) (24) 

into the net second - order contribution: 

E 2 = Eg N (n) + El°(n) + E™ (p) + E s h °{p) . (25) 

Here, (n) and (p) stand for neutronic and protonic contributions. See Eqs. (11)-(14) for 
the definitions of the different terms appearing in these two equations. This absorption is 
expected to have two major effects. Clearly, if such an absorption is possible, the factor a 
(see Eqs. (1) and (2) for definition), should be expressible as a function of neutron number, 
proton number, or some combinations thereof. Before discussing the possibility of absorbing 
fourth - order terms into second - order terms for a Woods - Saxon potential, we demonstrate 
the existence of such a functional form for the simple Harmonic Oscillator potential. 



A. The Harmonic Oscillator Potential 



The harmonic oscillator (HO) potential provides a unique opportunity to investigate the 
details of the WK expansions analytically. Therefore, first we consider the simplest form 
of the HO potential, without spin - orbit interaction. It can be shown that for the HO 



potential, the different WK corrections are given by 



10|: 



E < = ^sr (26) 
E > - m < 27 > 

where, A is the chemical potential, determined as described earlier, and u is the oscillator 
frequency. For the HO potential, assuming degeneracy of 2, the particle number (see Eq. 
(7)) is given by: 

1 / A \ 3 1 / A 



v =sU; -iU> (28) 
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This equation is cubic in X/(hu), and in principle can be solved exactly. Here, however, 
we take an alternative and physically more transparent approach, wherein, we express A as 



mm 



A = A + A 2 + A 4 (29) 

where, Xj is correct up to order HP. Starting from the Thomas Fermi expression for the 
chemical potential, and noticing that the normalisation is true order by order, we get the 
following expression for chemical potential, correct up to H 4 : 

A = |(3AT) 1/3 + i (3AT)^ 1/3 | ^o; (30) 

This, along with the second and fourth - order WK corrections to energy (see Eqs. (11), (12)), 
yields 

i 17 (M 2 ^ 

a wxfTK (31) 

where, X p and A n are chemical potentials for Z protons and N neutrons respectively. Further, 
notice that the neutron and proton numbers can be written as: 

N = ^-^A and Z = ^—^A (32) 

A = N + Z being the mass number of the nucleus, and I being asymmetry parameter, 
defined as I = (N - Z)/A. We obtain, 

a=i -^sr^ 2/3 ( i+ H < 33 > 

where, the terms up to the order A~ 2 ^ 3 are retained, and the expansion in I has been carried 
out only up to second - order in I. It can be therefore seen that the factor a can indeed 
be written as a function of mass number and J, implying that it is in principle possible, at 
least in the case of HO potential, to absorb the fourth - order WK corrections to the energy 
into the second - order WK corrections. 

To understand the behaviour of a with respect to /, we plot the factor a as a function 
of / in Fig. IHJ It is seen that the factor a has a very regular behaviour with respect to 
asymmetry. There are points stacked at a given value of /, with groups of points placed 
symmetrically with respect to them. This regularity persists over the entire range of / 
values. 
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FIG. 8: The factor a for Harmonic Oscillator potential, without spin - orbit interaction. Only a 
small portion of the asymmetry scale has been presented. 

B. Woods - Saxon potential 

Next, we investigate the factor a for the Woods - Saxon potential. In order to achieve 
this, we choose a set of 2171 known nuclei [25j with Z > 5. Spherical symmetry is assumed. 
The nuclear, spin - orbit and Coulomb potentials have been taken as defined in Ref. jsj. 
The full Wigner - Kirkwood calculations up to the fourth - order in h are carried out for 
these nuclei, and the exact values of the factor a are obtained. These are then plotted as a 
function of the asymmetry parameter / in Fig. |9j The figure exhibits that the factor a has 
a very regular behaviour as a function of asymmetry. In order to understand the detailed 
structure of the factor a, we plot the same results with a greater resolution in Fig. [101 

A remarkable and regular pattern emerges from the plots. In comparison with the case of 
the HO potential, the pattern is inverted. The pattern consists of 'fan like' structures. There 
are groups of points stacked exactly along vertical lines, as indicated in Fig. ITUl accompanied 
by symmetrically placed, slanting groups of points. All these groups of points constitute 
nearly perfect straight lines. This is in contrast with the case of HO potential, where the 
lines were curved. 

A closer examination of the behaviour of the factor a reveals several interesting features. 
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FIG. 9: Factor a as a function of asymmetry for a Woods-Saxon potential. 

To understand them better, we shall first enlist the nuclei appearing in a particular 'fan' 
structure. We shall designate the slanting lines appearing in the 'fan' structure as 'rays'. 
Thus, each fan structure has a number of rays in it, symmetrically placed with respect to 
the vertical line, defined by a particular ratio, (N — Z)/A. For example, let us consider 
(N - Z)/A = 1/11. This fan structure has 22 Ne, 33 P, 44 Ca, 176 Hg, ... along the vertical 
line. The first ray to the right of this line contains nuclei like 20 F, 31 Si, 42 K • • •. The second 
ray to the right of the vertical line consists of the nuclei like 40 Ar, 51 V, 62 Ni etc. The first ray 
to the left of the vertical line consists of 35 S, 46 Sc, 57 Fe etc. Whereas, the second ray to the 
left of the vertical line consists of 37 C1, 48 Ti, 59 Co etc. The heavier nuclei in this sequence 
are towards the bottom of the pattern. The value of a is therefore, inversely proportional 
to the mass number. Thus, it is expected that in the limit of A — > oo, the a values will 
approach some constant value, say, a , which is approximately 1.125, according to the figure 
above. 

Considering these observations, we propose the following parametrisation for the factor 



a: 



a x N - Z 
a = a ° + ~A + a<1 — A — 



a 3 



N-Z 
A 



(34) 



where, a/s are adjustable parameters. Considering all the 2171 nuclei (see above), we carry 
out a least squares fit to determine these parameters. The fit turns out to be exception- 
ally good, with rms deviation 1.09 x 1CT 3 . The values of the parameters are: a —l. 12761; 
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FIG. 10: Fate a as a function of asymmetry. 

a i=2. 26744; a2=-0. 02659 and ^3=0. 29987. The difference between the exact and the corre- 
sponding fitted a values is plotted in Fig. [HJ indicating that the agreement is almost perfect, 
and that the phenomeno logical formula that has been proposed here is indeed robust, for 
all the mass regions. 

We shall now investigate the deformation effects particularly with reference to the factor 
a. In order to achieve that, we once again consider the set of 561 even - even nuclei (see 
Section 3), with deformation parameters obtained as described before. The calculation of 
binding energies requires the shell corrections, pairing energies and the liquid drop energies. 
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FIG. 11: The difference between the fitted and the corresponding exact values of a. 



The shell corrections require averaged energies, which are calculated here using the WK 
expansion. Here, we consider the WK expansion only up to second - order, and simulate the 
effects of fourth - order through the factor a (Eq. ( 134")) ). This defines the averaged energies 
and hence the shell corrections completely. The difference between the shell corrections 
thus obtained and the 'exact' shell corrections is found to be indeed small, the maximum 
deviation being of the order 150 keV, implying that the factor a obtained merely by using 
the spherical nuclei works very well for deformed systems as well (with both deformation 
parameters and f3^). This observation is indeed of great practical importance. 

With these approximate shell corrections, we make a re-fit to the liquid drop parameters. 
Comparison between the liquid drop parameters as reported in Section 3 and the ones ob- 
tained with the approximate shell corrections is presented in Table 3. It is indeed gratifying 
to note that the liquid drop parameters obtained in the two cases are almost identical, and 
so is the rms deviation of the calculated binding energies with respect to experiment 25 ]. 
This substantiates the validity of the parametrisation of a. 

To test the robustness of the parametrisation of a further, we calculate the constants <x,'s 
in Eq. (El using just four nuclei ( 40 Ca, 100 Sn, 146 Gd and 208 Pb) instead of 2171 nuclei as 
described above. It is found that the numerical values of the constants practically remain 
the same. To test the validity of these parameters, the liquid drop parameters are re-worked 
employing the new values of a/s. It is found that the liquid drop parameters thus obtained 
are practically equal to the ones reported in the right most column of Table 3. 
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TABLE III: Values of the liquid drop parameters obtained through the x 2 minimisation for 'exact' 
and 'approximate' shell corrections denoted by 'E' and 'A' respectively. 



Quantity 


'E' 


'A' 


a v 


-15.435 


-15.421 


k v 


-1.875 


-1.873 


a s 


16.673 


16.580 


k s 


-2.430 


-2.432 


O'CUT 


3.161 


3.295 




1.219 


1.221 


c 4 


0.963 


0.953 


Wl 


-2.763 


-2.652 




3.725 


3.659 


rms 


0.610 


0.607 



We close this section, by concluding that the absorption of fourth - order Wigner - 
Kirkwood corrections into the second contributions is reliable, and can be used in large scale 
mic-mac calculations. The absorption also has the advantage of reducing the numerical noise 
that might arise in the higher order derivatives of the potentials. 

VI. SUMMARY AND CONCLUSIONS 

The semiclassical Wigner - Kirkwood h expansion of the one - body partition function has 
been employed instead of the Strutinsky averaging scheme to calculate the shell corrections 
within the framework of a mic-mac model. The microscopic part of the energy also contains 
pairing contributions that are obtained using the Lipkin - Nogami scheme. We have improved 
the macroscopic part of the model as compared with the one used in our previous work Q] 
by including the curvature term as well as the Wigner contribution. With just ten adjustable 
parameters, our model reproduces the binding energies of 561 even - even spherical and 
deformed nuclei with rms deviation of 610 keV. We have tested this new mic-mac model 
near the proton and neutron drip lines as well as in the superheavy region. Our present 
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calculations show that the mic-mac model proposed in this paper reproduces remarkably 
well the recent experimental results in these exotic scenarios. 

Further, a systematic study of the ratio of the fourth - order and second - order Wigner - 
Kirkwood energies has been carried out. We find that the ratio of these two energies behaves 
in a very systematic manner. We have shown that this ratio can be parametrised accurately 
by a simple expression, implying that the fourth - order corrections can be absorbed into the 
second - order contributions in a very simple way. We have checked that using this simple 
procedure, we recover practically the same parameters of the macroscopic part, without 
deteriorating the quality of agreement achieved with the full Wigner Kirkwood calculation 
including explicitly the fourth - order contributions. Therefore, this simplified calculation 
of shell corrections can be used confidently in the large scale mic-mac calculations that we 
plan to carry out as the next step. 

Finally, we point out that there is still some room for improving our model particularly 
in two specific directions. On the one hand, the full blocking procedure in the pairing 
calculations of odd - odd, odd - even and even - odd nuclei, that may be particularly 
relevant for spherical nuclei, has to be introduced. On the other hand, refinements in the 
mean field Woods - Saxon potential and in the distance function are still needed to study 
with our model not only neutron rich nuclei, but also fission barriers. This would require 
large scale calculations with the model, for which, the simplification proposed above may be 
very useful. 
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